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Abstract 

We study the evolution of flavor-singlet, light-cone amplitudes in the soft-collinear effective 
theory (SCET), and reproduce results previously obtained by a different approach. We apply our 
calculation to the color-singlet contribution to the photon endpoint in radiative T decay. In a 
previous paper, we studied the color-singlet contributions to the endpoint, but neglected operator 
mixing, arguing that it should be a numerically small effect. Nevertheless the mixing needs to be 
included in a consistent calculation, and we do just that in this work. We find that the effects of 
mixing are indeed numerically small. This result combined with previous work on the color-octet 
contribution and the photon fragmentation contribution provides a consistent theoretical treatment 
of the photon spectrum in T — > ^X. 
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I. INTRODUCTION 



The soft-collinear effective theory (SCET) 0,0,0,0] is a systematic treatment of the high 
energy hmit of QCD in the framework of effective field theory. Prior to the introduction of 
SCET this limit of QCD was subject to intense study using various other approaches includ- 
ing all order perturbative methods . Some of these classic calculations have been revisited 
in SCET and their results reproduced 0,0. The effective theory approach, however, goes 
beyond the approximations upon which many of the previous calculations rely. In particular 
it is straight forward to include power corrections to any process, as was demonstrated in 
the context of color-suppressed B meson decays 0, which receive their first contribution 
at subleading order. In addition SCET naturally includes nonperturbative effects in the 
form of matrix elements of operators. This for example, gives a consistent explanation for 
the origin of the shape function in semi-inclusive B meson decay. In this article we study 
the radiative decay of the Upsilon, and revisit another classic result: namely the evolution 
equation for light-cone wave functions, also known as the Brodsky-Lepage equation. 

At first sight it may seem strange to be discussing a heavy quarkonium system in the 
context of a high-energy effective theory. It is, however, the final state of radiative Upsilon 
decay which can, in a certain region of phase space, be described by SCET. To describe the 
T system, which is a boundstate of a heavy b quark and b quark, we need to consider a 
different limit of QCD: the non-relativistic limit. This is sensible because the large mass of 
the b quark ensures that the the typical relative velocity v of the b and b in the T is small 
V ~ 0.3 allowing for a non-relativistic expansion. Furthermore the production and decay a 
bb pair can be calculated perturbatively. In the earliest works on quarkonium, the f — > 
limit was always taken, allowing all the non-perturbative dynamics to be isolated into the 
wavefunction at the origin. This approach is now called the color-singlet model, since in an 
effective theory picture it corresponds to keeping only those operators that create/annihilate 
the bb in a color- sin glet configuration. With the advent of a non-relativistic effective theory of 
QCD (NRQCD) Sll^, this simple picture is replaced by a systematic expansion in operators 
that scale as higher and higher powers of v, where some of the time the quarkonium state 
can be produced/annihilated in a color-octet configuration. 

The theoretical picture of radiative Upsilon decay that emerges from these considerations 
is quite rich. Over some of phase space the decay is described by the annihilation of a bb pair 
in a color-singlet ^5*1 configuration into a photon and a pair of gluons with invariant mass 
on the order of the Upsilon mass. This is well described by an operator product expansion 
based on NRQCD. However, the situation is more complicated as the photon energy reaches 
its maximum. In this region of phase space the pair of gluons form a coUinear jet back-to- 
back with the photon, and there arises a possibly large contribution from the annihilation 
of the bb in a color-octet configuration into a photon back-to-back with a single gluon. Since 
the decay products in this "endpoint" region are jet-like (i.e. their energy is large relative to 
their invariant mass) the appropriate effective theory to describe the dynamics of the decay 
products is SCET. The T is still described by NRQCD. 

The radiative decay of the Upsilon is of particular interest since it allows for a measure- 
ment of the strong couphng constant ll[ 13,113, 13 • Furthermore the differential decay 



rate as a function of the energy fraction z = 2E^/Mr has been measured, and in each case 
found to be softer than the QCD predictions. In a series of recent papers [l^, [13, [l3| this 
decay has been studied using SCET to describe the endpoint region of the decay rate. First 
in Ref. jl3], the large Sudakov logarithms for the color-octet channels were resummed for 
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the first time using SCET. In subsequent papers we analyzed the color- sing et decay in 
the endpoint region. This was calculated by Photiadis However, in Ref. Il6| we ignored 
the possibility of a jet of a light quark and anti-quark in the final state, and we reproduced 
Photiadis' results in this limit. The qq final state has a zero tree- level matching coefficient 
in the effective theory for this process, but it can be generated by mixing with the gluon jet, 
and so it must be included in a consistent calculation. 

The main result of this work is the derivation within the SCET framework of the evolution 
equation for matrix elements of coUinear operators that describe the gluon and quark jet final 
states in radiative Upsilon decay near the endpoint. As a consequence of the factorization 
of soft physics from collinear physics the evolution of these matrix elements of helicity- 
zero, flavor- and color-singlet collinear operators is quite general, and should hold for any 
collinear final state produced from the vacuum. This was pointed out by Photiadis. Thus, 
the evolution equation for the matrix elements we are concerned with should be the similar to 
that of the pion lightcone wave function which was first considered in Ref. |l9ll2o|. However, 
the full mixing is not incorporated in those works, since the pion is a flavor non-singlet. The 



flavor singlet case was done by Chase [21| in the context of quark- ant iquark and gluon-gluon 
jet production in photon-photon collisions. We reproduce those results using SCET. In 
Sec. II we quickly review soft-collinear effective theory, in Sec. HI we introduce the collinear 
operators that arise in radiative Upsilon decay, in Sec. IV we calculate the renormalization 
group equation that governs the running of these operators, in Sec. V we use the results 
of the previous section to give the resummed rate for radiative T decay, and in Sec. VI we 
conclude. 



II. REVIEW OF SCET 

We begin with a short review of the parts of SCET that are relevant to this calculation. 
In particular, we are only concerned with SCETj ji^l, which describes the interactions of 
collinear and ultra-soft (usoft) degrees of freedom. In this theory collinear particles have 
momenta whose lightcone components scale as p = ~ 1, A), where Q is a 

large energy scale, and A ^ 1 is a small expansion parameter. In SCETj A ~ a/Aqcd/Q 
so that the typical invariant mass of an SCETi collinear particle is ~ Q^q,ct>- Usoft 
particles have momenta which scale as A; = (A;"*", fc^) ~ Q(A^, A^, A^), so that the typical 
invariant mass of a usoft particle is fc^ ~ Aqcd^- The usoft degrees of freedom interact with 
the collinear particles without taking the collinear particles off-shell by more than ~ QAqcd- 
Furthermore it is only the plus component of the collinear momentum that a usoft particle 
can change. 

Here we are interested in the differential decay rate for T — > 7 + X as a function of the 
photon energy restricted to the region where 2E^ ~ My — (9(Aqcd)- In this regime the final 
state hadrons have a lightcone momentum componenet of order Mx, and invariant mass 
of order a/Mx(Mx - 2E^). Clearly the jet can be described with SCETj, where Q = Mr 
and the power-counting parameter is A = ^Jl — 2E-y/Mf. The T particle can be treated 
m NRQCD where large fluctuation about the heavy quark mass are integrated out, 

leaving only modes with momentum of order mv"^ ~ Aqcdj where v ~ 0.3. These usoft 
modes can interact with both the heavy quarks in the inital state and the collinear particles 
in final state. 

By matching QCD onto SCET the large scale Q is integrated out. In practice, the 
matching procedure is to calculate matrix elements in QCD, expand them in powers of A, 
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and match onto products of Wilson coefficients and operators in SCET. Thus it is important 
to be able to deduce the SCET operators which can arise at a given order in A. Field theory 
generally allows all operators that are consistent with the symmetries of the theory. As 
explained in detail in Ref. j^, the symmetry of SCET which restricts the operators that 
can arise is gauge invariance. Specifically, SCET is invariant under two types of gauge 
transformations: collinear and usoft. Under coUinear gauge transformations the usoft fields 
remain invariant, while the collinear fields transform in the usual manner. Under usoft 
gauge transformations the usoft fields transform in the usual manner, and the collinear 
fields undergo a global color rotation. 

The collinear fields in SCET are the fermion field ^n,p, and the gluon field A!^^. These 
fields are labeled by the lightcone direction n^, and the large components of the lightcone 
momentum (n ■ q,q±)- The femion field contains a term C,np that annihilates particles, and 
a term ^n-p that creates antiparticles. For the construction of gauge invariant operators we 
will find it convenient to make use of the SCET collinear Wilson line, 



exp (^-gs^n ■ An,g{x)^ 



'- perms 



(1) 



The operator projects out the momentum label [3| of fields that sit to the right of the 
operator. We will use the convention that only acts on those fields in the square brackets. 
Generally [/ (P) J, ■ ■ ■ J„ 0pi ■ ■ ■ 0p„ ] = / (n ■ pi + ■ ■ ■ + n ■ p„ - n ■ gi ^ ■ ^m) ■ ■ • , 

where P = n-V. The conjugate operator acts on fields that sit to the left of the operator, 
and projects out the sum of labels on conjugate fields minus the sum of labels on fields. In 
the usoft sector there is a usoft fermion field qus, and a usoft gluon field A!^^. 

Operators in SCET are constructed such that they are gauge invariant under both 
collinear and usoft gauge transformations. For example, under collinear-gauge transfor- 
mations ^n,p Unin,p and Wn UnWn, SO the Combination 

WlUp (2) 

is collinear-gauge invariant. Furthermore it is convenient to introduce a delta function which 
fixes the labels of the combination of fields above: 

Xn,Lo = [Su;,pW^^n,p\, (3) 

where it is understood that we will include a sum over uj for each Xn,uj in an operator. The 
Wilson coefficient will in general depend on the label momentum uj, which will result in 
a convolution of the short distance coefficient with the SCET operator. The combination 
above still transforms under a usoft-gauge transformation Xn,uj — ^ V{x)xn,uj- 
The collinear-gauge invariant field strength is 



7 

/^/l!/ 

n 

9s 



(4) 



where 



= -V + Vl + -m ■ D, (5) 

and iD'^ = id^ + gsA^ is the usoft covariant derivative. Note that Gi^ is not homogeneous in 
the power counting. The leading piece scales like A, and is given by \PB^ = fiyG^, where 
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the perp subscript on B indicates that the /i index only has support over perpendicular 
components. Simplifying and including a label fixing delta function, we obtain 



5^ = — 



(6) 



Under usoft gauge transformations V{x)B'^^V''{x). We use these objects to build 

the operators we need to match onto SCET at the endpoint of the T — X7 spectrum. For 
further examples the reader is referred to Ref. 



III. SCET OPERATORS 

We begin by matching the QCD final states onto SCET operators. Since we are interested 
in the color- and fiavor- singlet, helicity-zero operators, we have the SCET operators 

Oq{uJl,UJ2) = Xn,u^^2^n,c.^, (7) 

where g'^^ = g°'^ — {n°'n^ + n^n"')/2. We have introduced an additional factor of V into 
the gluon operator so that both of the above operators have the same energy dimension. In 
addition, both operators scale as in the SCET power counting. They are the complete 
set of leading color-singlet operators. Each of the operators are convoluted with a short 
distance coefficient Tg/g{uJi,uj2), which is determined by matching onto QCD. 

Matrix elements of the operators in Eq. ((7j) are non-perturbative functions of the labels 
Ui and UJ2- Consider the matrix element of a coUinear final state and collinear initial 

Xn,w\ 2 Xn,ui2 \-^n,p' )• (8) 

This can be simplified by introducing uj± = u;i ± 0^2 and p± = p ± p', and using momentum 
conservation 

<Pri{x+) , (9) 

where is the light-cone amplitude (LCA) for the transition X ^ T ^ and is by definition 
dimensionless. This last requirement on 0;ri forces us to introduce the constant /Cjfi which 
is process dependent and possibly dimensionful. In Upsilon decay Kj = M^. To arrive at 
the last line of the above equation we extend the sum over discrete to an integral over 
continuous UJ+ and define a:+ = a;+ /n ■ p_ . As a result all sums over u;+ are converted to 
integrals over x. In Appendix ^ we show how this is done using type (a) RPI as defined in 
Refs. 113. 

Two specific choices of initial and final state are familiar. If we choose the incoming 
and outgoing state to be a proton with momentum p, then is related to the parton 
distribution functions. If, however, the incoming state is the vacuum, and the outgoing 
state is a meson with momentum p, then is related to the light-cone wave function of 
the meson. 
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In the case of T ^ 7 + X in the large photon energy regime the QCD amphtude for 
bb{l,^Si) — > 700 matches onto a convolution of a short-distance Wilson coefficient and an 
SCET current 

where 

Jii,^s,M^f^) = xlpA-^Vr[Bl6^,T.^B^] . (11) 

The NRQCD fields ipp and X-p annihilate the heavy quark and antiquark fields, respectively. 
From now on we will drop the (Ij^Si) label. Note we correct a typo in Ref. [l6| where the 
Kronecker delta has the incorrect label operator. The Upsilon contains no collinear quanta, 
so the current factors into an usoft piece containing the heavy quark spinors and a collinear 
piece containing the trace over the SCET gluon fields. The usoft fields can not "talk" to the 
collinear fields due to color-transparency. We have simplified the above expression by fixing 
the momenta to be those of the the particular decay we are interested in. Strictly speaking 
this can only be done after taking the matrix element of the operator above between external 
states, which is given by 

=^ {X^\xlpA-(T^iJp\T) j^^dxTg{x;fi)(Pg{x;fi). (12) 

The outgoing state, Xc, is a jet with total momentum p~ = Mr, fixed by the mass of the 
decaying Upsilon. The kinematics are similar to the meson light-cone wave function, and 
we therefore expect the running of the collinear matrix element (f)g{uj] fi) which appears here 
in Upsilon decay to be the same as the running of the light-cone wave-function of a meson. 
Note the usoft matrix element does not run 12511. 



IV. RUNNING OF OPERATORS 

In SCET large logarithms are summed using the renormalization group equations (RGEs). 
In the case we are interested in there are two LCAs, the matrix elements of the operators 
in Eq. ((7j), and they mix with each other. This will result in a coupled differential equation. 
In addition the LCAs under consideration is functions of the momentum fraction x, which 
makes the RGE an integro-differential equation. 

The bare SCET operators, denoted by a zero superscript, are related to the renormalized 
operators through a counterterm: 

0^a\^) = J dyZab{x,y;fi)Ob{y;fi) a,b = g,q, (13) 

where the bare operator does not depend on the scale /i. Differentiating this equation with 
respect to and using the identity 

J dz Zab{z,x; iJ,)Z~^^{z,y; n) = SbcS{x - y) . (14) 
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we obtain 



- / dy [ fi-^Zab{x,y]fi) ]Ob{y]fi) 



dyich{z,y]^)Ob{y]^) , (15) 



d 

dy Zab{x, y; /i) fi—Ob{y; fi) 
djjL 

=^ /^^C'c(z; /i) 



with 'jcbiz,y; n) the anomalous dimension. 

The running of the short-distance coefficient is obtained as a consequence of the scale 
independence of the full theory current. For example differentiating both sides of Eq. 
with respect to yU gives zero on the left-hand side, which then gives a relationship between 
the running of the operators and the coefficient function. Generally the QCD current is 
matched onto the full set of operators given in Eq. (|7j), and differentiating we obtain 



d f 

= /i— J dxTaix; fi)Oaix; iJ,) 



(16) 



dx 



dx 



Oaix; n)iJ-^Taix; fx) + Taix; fj)fi-^Oa{x\ /i) 
djjL djjL 

d f 
Oa{x] /i) fi—Ta{x; f^) - J dy Ta{x] fx)-fab{x, y; fx)Ob{y] /x) 



dx Oa{x] /i) 



/i|^r.(a;;/.) 



dyTb{y]n)-fba{y,x;fi) 



0, 



where we use the result of Eq. fll5|) in obtaining the penultimate expression above, and inter- 
changed the X and y (and a and b labels) in the second term to obtain the final expression. 
Since this must hold for any value of x, this equation implies an RGE for the coefficient 
function 

d f 

/^^ra(x; Ai) = y dy Tb{y; fi^baiy, x; /i) . (17) 

The renormalization Z can be calculated in perturbation theory. In dimensional regular- 
ization [MS scheme) we obtain to 0{as)'- 



Zab{x, y; fx) = 6abS{x - y) + - ^''}^^ Pab{x, y) . 



e 271 



In order to satisfy Eq. p4|) we must have 



Zabi.^^ /^) = ^abK^ -y)- -'^^^Pab{y, X) , 



(18) 



(19) 



where on the right hand side of the equation above the x and y have been reversed. We now 
proceed to calculate Z for the matrix element which arises in Upsilon decay. 

The Feynman diagrams which are needed to calculate Z are shown in Figure ^ while 
the Feynman rules for the operator vertices are given in Appendix^ We show only those 
diagrams that are non-zero in dimensional regularization where the infrared is regulated by 
choosing pi = pi = and Pi-p2^ 0, with pi 2 the momenta of the final state particles. The 
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FIG. 1: One loop renormalization: a) glue to glue, b) quark to glue, c) glue to quark, d) quark to 
quark. The quark and gluon lines all represent collinear particles. 



divergent piece of the amplitude for each set is 

dy Aab{x, y)Ta{x; fi)(j)b{y; /i) 





27re 


^gg{x,y) 


= -Ca 


^gq{x,y) 




^ggi^^y) 


2 I 


\q{x,y) 





+ y'^ 



1 

+ - 



1 



L(l + x)(l-y) 2 ix-y)+\ 
l-x + 2y 



hJ = q,9 , 
eix-y) + 



(20) 



y^-y 



6{x-y) 



(l-x)(l + y) 
1 — X 



Q{y-x) - 



y^-y 



l-y + 2x)Q{x -y) - 



y^-y [ ' 



'1 


— X 






1 






x-y J 



Q{x-y) + 



y-^-y 



In obtaining these expression we made use of the property Ta{—x) = Ta{x), which is a 
consequence of the invariance of the product of operator and coefficient function under 
charge conjugation. These divergent amplitudes are canceled by the renormalization Zat- 
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First we invert Eq. ()13|) and take the matrix element 

'dy{Ol'\y))Z-,\y,x) 



(21) 



where Zf is the renormahzation factor for the fields in the operator Oa- Expanding this 
to first order in dimensional regularization where Zf = 1 + as(/i)5//(27re), and using the 
expression in Eq. (fT^ we obtain the equation which fixes the Pab- 



Pab{x,y) = Aab{x,y) + 6f6abS{x - y) . 
From this we calculate the one loop expression for the anomalous dimension 

as{jj) 



^ab{,x,y]ii) 



IT 



-Pabix, y) 



(22) 



(23) 



and substitute it into Eq. (fTTj) to obtain the one loop RGE 



H-^Taix;fi) 
an 



7T 



dy Tbiy;n)Pbaiy,x) 



dy Tb{y;fi) 



Aba{y, x) + 6f6abS{x - y) 



(24) 



With this result in hand we can solve the RGE by diagonalizing. The first step is to 
expand the coefficient funtions in a basis which is diagonal under the convolution with the 
Pab- This basis is provided by the Gegenbauer polynomials 



r,(x,/i) = 5^4'^)(/i)(l-a:^)C^/_^,(x) 



(25) 



n odd 



n odd 



where the restriction to odd n is required by Bose symmetry for the gluons. Substituting 
these expansions into the evolution equations yields coupled ordinary differential equations 



where 



(n) 
Iqq 



'gq 



'qg 



'gg 



d 

djj \ a 



Cf 
1 



(n) 
ttq 

in) 

g 



TT 



n) (n) 

Tqg Tgg 



n+l 



(n) 
ttq 

(n) 



_(n+l)(n + 2) 2 



9 ^ 7 



(26) 



(27) 



n 



3n + A 



3 ^(n + l)(n + 2)' 
n"^ + 3n + A 
n{n + l){n + 2){n + 3) 



Ca 



+ 



n{n + l) {n + 2){n + 3) 



1 
6 



n+l 



^ 1 



i=2 



-rii 
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The RGE in Eq. fl26p can be diagonalized through a simlarity transformation resulting in 



/i 



where the matrix A is diagonal and has eigenvalues 



^W? + 7j?±A], with A = ^ 



(28) 



(29) 



The eigenvector a*^"-* is 



An) 



An) 



(n) (n) (n) / > (n) (n) 

Ctg 799 - «9 \^+ 



(30) 



The diagonalized RGE is simple to solve, giving 

a,(M)T"'^±/^° 



ai")(/^) 



as(/i) 



ai^^M) 



where /3o = 11 — 2n//3. The equations above can be inverted to obtain 



(31) 



(32) 



\{n) (n) (n) (n) 



99 



99 



We can now include the running of the coefficients to get a result for the resummed gluon 
and quark coefficient: 



as(M) 



1 -2aV"VA) 



99 



_ la(")(M) 
A ^ ^ 

-2a(."V/3o 



as(M) 



(33) 



A7,^") 



99 



So far our results have been general, and can be used for not only Upsilon decay, but any 
process with helicity-zero, flavor- and color-singlet wavefunctions. The process dependence 
will come in the boundary conditions. For Upsilon decay, the matching coefficient for the 
quark operator is zero at leading order, while the matching coefficient for the gluon operator 
is a constant k. We will normalize our matrix element so that k = 1. Having expanded the 
matching coefficients in Gegenbauer polynomials we determine 



a(")(M) 
4")(M) 







(34) 



where 



An) ^ l)(n + 2)(n + 3) 

•^5/2 9(n + 3/2) 



(35) 
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is the normalization of C^^i{x). 

Using the relations in Eq. ()32|) we determine the initial conditions for the componenets 
of a: 



These can be substituted into Eq. to obtain the final result: 



A 



(n) 

7+ 



aW(M), 



-7- 



(n) 



a(")(M) 



where 



(n) 



An) 



A 



(36) 



(37) 

(38) 



(39) 



V. RESUMMED RATE 

The decay rate is proportional to the imaginary part of the forward scattering amplitude 
T. The expression for this ampitude was derived and given in Eq. (59) of Ref. [l^^, 



ImT(2;) 



dx 



2M- 



T 



M2 



H{x) / dtS{t, ^)ImJ[Mx, t + M(l - 2); /x] , 



(40) 



where if is a hard coefficient, S{1'^) is the color-singlet shape function j2 

dx 



s{i- 



■,t 



■e-^'"^"/'(T|T[^V«X-p](a:-)[xLp,a.^p'](0)|T) 



(41) 



and lm.J[Mx{i^ + M(l — 2)); /i] is the imaginary part of the jet function, 
(OlTTr [Bl5{u - m-I?+)5^] {y) [B^diJ - in-V^)B^] |0) 



2i5{uo - uj') 



d^k 



-'''■yj{uj,k+;fi) 



(42) 



where the labels iv and co' are continuous and in-V = V + in-d as discussed in Appendix 1X1 
Since Ref. did not consider mixing, this was the only jet function. We now generalize 
this to 



{Q\T[Oa{uj; y)0,{uj'- 0)]|0) = 2iM 6{uj - uj') 



d^k 



-*V,,(u;,A;+;/i) 



(43) 



^ Here we fix a typo in that equation. 
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FIG. 2: Feynman diagram for the leading order gluon jet function. 




FIG. 3: Feynman diagram for the leading order quark jet function. 



where a,b = g, q. The hard coefficient gets modified to be 



(44) 



If a 7^ 6, we get no contribution at this order in perturbation theory. When a = b = g, we 
get exactly what was considered in Ref. pictured in Fig. |21 The imaginary part of the 
jet function in this case is 

lmJgg{uj,k+;fi) = ^e{k+). (45) 



We now need the imaginary part of the quark jet function, picture in Fig. El The result is 

N, M^-uj'^ 1 



ImJgg(a;,/c+;/i) 



m - 1 M2 Stt 



(46) 



Expanding the matching coefficient in Gegenbauer polynomials the integral over u = Mx 
in Eq. (^0)) may be be carried out, giving 



1 „ J~l 



X 



2\2 



n odd 

16 



n odd J 5/2 



for the gluon jet function, and 



[\x{l~x^)Tl{x,fi) = V K"^/")]' f dx{l-x'')Cl'\x)Cl'\x) 

n odd "^-1 



16 



f(") (n)2 

9 ^ r.Hp 

n odd U5/2J 



"q" 



^(^)<V/3o_^(^)2A-//3o 



1 2 



(47) 



(48) 
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for the quark jet function, where we have defined 

asifi) 



r(/i) 



and 



asiM) ' 
n + l)(n + 2) 



An) 

•^3/2 n + 3/2 



(49) 



(50) 



is the normahzation of Cn^^(x). Using the results of Ref. [l6j, the differential decay rate is 



1 dTr 



Tq dz 



e{Mr-Mz)— ^ 



9 ] An) 



7? V(/X,)2A-7/3o _ ^i"V(/X,)2A-7/3o 



1 2 



(51) 



n odd 



qf(") (n)2 



where /Xc = M\/l — 2: is the coUinear scale. This results differs from the result in Ref. [fi 
We agree with Ref. [l^ up to Eq. (5) of that paper, up to some typos. However, following 
the method outlined in that paper, we still arrive at the first line of Eq. (|5ip. Both Eq. (|5ip 
and the result in Ref. jl8|] reduce to the rate calculated in Ref. when the mixing is 
turned off. The second line of Eq. (j5H) comes from the quark jet function. This adds a small 
contribution to the total resummed rate. 

While the logarithms that are summed in Eq. ()5ip are important at large this formula 
should not be trusted away from the endpoint. In order to match our resummed result onto 
the leading order result, we will interpolate between the two using the formula. 



1 dV 



int 



Fn dz 



where l2 



1 dV 



LO direct 



2-z z{l 



dz 



(2 



1 dVth 
Fn dz 



z) 1 
-t4 + 2- 



+ 



1 dVr 



Fn dz 



ln(l 



and 



Fn= 



32 2 2 (^l^-p^^^-p^-P'^^^p'l^) 



(52) 



(53) 



(54) 



The term in brackets in Eq. ()52|) vanishes as ^ — 1, leaving only the resummed contribution 
in that region. Away from the endpoint the resummed contribution combines with the —z 
to give higher order in as{M) corrections. This is clear from Eq. (jKT|l . There are important 
corrections to this result due to fragmentation at low z |Q| . However, since we are interested 
in the large z region, we will neglect them in the following. There may also be large color- 
octet corrections to the rate in the endpoint region 0, 17, |2^, which we will also neglect 
for the now. We also compare our result to the resummed result where the mixing has been 
turned off, using 



1 dVr 



Fn 



dz 



e{Mr-Mz)— J2 



(n) 

n odd J 5/2 



Qs(/ie) 

a.(M) 



(55) 
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FIG. 4: The color-singlet rate. The dotted curve is the tree-level direct rate. The solid curve is the 
interpolated resummed direct rate. The dashed curve is the resummed rate with the mixing turned 
off. 

in place of the full resummed result in Eq. (jH21)- 

In Fig. |3]we show the color- singlet, interpolated resummed rate, Eq. (|52j) show as the 
solid line, compared to the leading tree level color-singlet result, Eq. (|^. shown as the 
dotted line. As can be seen, the resummed rate turns over and decreases near the endpoint. 
Also shown in Fig. |^ is the interpolated resummed rate with the mixing turned off, Eq. (|55|) . 
This is the same as the results of Ref. The result without mixing is a fairly good 

approximation to the full result. We also show the contribution coming from the quark jet 
alone as the dot-dashed line. 

VI. CONCLUSIONS 

Radiative Upsilon decay at maximal photon energy is characterized by a photon recoiling 
against a jet of coUinear particles. Thus SCET is the appropriate effective field theory to 
study this kinematic situation. The lowest order color-singlet QCD diagram for this process 
has the Upsilon decaying to a photon and a pair of gluons. In a previous pair of papers 
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[l6l |. we used SCET to investigate the endpoint behavior, summing kinematic logarithms. 
However, we neglected the possible mixing of the gluon pair with a quark-ant iquark pair. 
The full calculation, including the operator mixing, had been presented in the literature by 
Photiadis As pointed out in Ref. [l^, the radiative Upsilon decay at the endpoint has 
the same evolution equations as the flavor- singlet light-cone wavefunction evolution. 

Therefore, we have calculated the flavor- and color-singlet, helicity-zero light-cone ampli- 
tude evolution using SCET, with the goal in mind of studying the photon endpoint spectrum 
in radiative Upsilon decay. We find that SCET does reproduce the evolution equations for 
the light-cone amplitudes presented previously in the literature 0|. When applying this 
to Upsilon decay, we however disagree with Ref. although numerically the results are 
similar. 

With the inclusion of the operator mixing, we have a complete, leading logarithm result 
for the color-singlet contribution to radiative Upsilon decay at the endpoint. Combi ning 
this with the leading logarithm result for the color-octet contribution at the endoint 3, ITf - 
and the photon fragmentation results at low z psi . 2^, we can hope to obtain an accurate 
prediction for the photon spectrum over the full kinematic range. 
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APPENDIX A: OPERATORS OF CONTINUOUS LABELS 

In this Appendix we explain the relationship of SCET operators defined using a discrete 
label to those defined using a continuous label. As a concrete example we consider the 
current in Eq. ()12|). which involves the gluon operator. The matrix element of the collinear 
operator in the first line is 

(Xe|Tr[5^5^,P^5^](x)|0), (Al) 

where we have made explicit the space-time dependence. The expression above is defined 
for a discrete label a;. However, we could write down an operator involving a continuous 
label 

(Xe|Tr[5^5(^, - in-V^)Bi\ (x)|0) , (A2) 

where in-V^ = V+ — in - d + in- d . Note that the sum over u in Eq. (fT^ is now replaced 
with an integral over Uc. The delta function must be understood as 

6{lJc — in-V^) = S^p^6{k — in-d+) (A3) 

where ujc = uj + k with ~ M discrete, and k ~ Aqcd continuous. The integral over ujc must 
then be understood as a sum over uj and an integral over k. The expression in Eq. ()A2|) can 
be expanded in powers of in-d/V ~ Aqcd/^, where the leading term is just the operator 
in Eq. ()A1|) . Thus the continuous operator is just the discrete operator plus high-order 
corrections. However, in an EFT it is only sensible to include higher order corrections in a 
leading order operator if all of the subleading terms run the same way as the leading term 
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(i.e., they all have the same anomalous dimension). This is only true if there is a symmetry 
which enforces this condition. In this case the symmetry is a specific reparameterization 
invariance known as RPI (a) |23l . 2^ . In essence this RPI is the statement that there is no 



unique way to decompose the large label momentum and the continuous residual momentum. 
This implies that reparamterization invariant operators must be built out of in-V, and that 
such an operator runs in a specific way. As a result any subleading operators that are due 
to an expansion of in-V in powers of in-d/V must have the same running. 



APPENDIX B: FEYNMAN RULES 

In this Appendix we give the Feynman rules derived from the color-singlet operators given 
in Eq. (|7j), which we repeat here: 

The fields B"^^ and Xn,uj are built using the coUinear Wilson line in order to obtain gauge 
invariant objects. We thus have an infinite number of Feynman rules encoded in each 
operator of Eq. ()B1|) . Here, we will give the corresponding Feynman rules necessary for 
calculating the anomalous dimension of the operators at one loop, namely the operators to 
order and gl- We will always define our momentum to be incoming. The Feynman rules 
are shown in Fig. El 

We begin with the gluon operator. We have 
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where uj± = uji±uj2, V± =V±V\ and the factor of 2 will cancel the Jacobian from changing 
from uJi^2 to uj±. The delta function 6{uj- — V-.) will constrain the sum of the momenta to 
be the total energy of the jet, which in our case is Mr- We are therefore only interested in 
the Feynman rule for 

0,(^_) = PTt[BI 5K - P+) Bi]g^^ . (B3) 
Expanding out the B'^ to leading order in Qs, we get the order Feynman rule 

iAf = ~M6^^ [6{uj+ + n-g,-n-g2) + 6{uj+ -n-g^ + n- g^)] 



^9^u 9± - ] [91 - • (B4) 



At order g] we get 
2 



Og(cJl,CJ2) = 2^n,pi-K^+-'P+)8{uJ_-V-)^n,p2 



^■9i J \\ n-{g2 + ga) J n-g2 \ n-{g2 + gs) J n-g^ } 
+ [(l,a)-.(2,/5)^(3,7)^(l,«)] 

+ [(1, a) (3, 7) - (2, /3) ^ (1, a)] . (B5) 
We can similarly rewrite our quark operator as 



=^ 0,icu_) = - V+) • (B6) 

This gives the order g^ Feynman rule 

tAf = i^n^^J{uj+ + n ■ gi - n ■ ga) , (B7) 
where again, the momentum is defined to be incoming. The order gl Feynman rule is 

MJ^^ = ig.s-^L^T'^in [S{uj+ + n ■ qi - n ■ q2 - n ■ g) - S{uj+ + n ■ qi - n ■ q2 + n ■ g)] . 
^ n ■ g 2 

(B8) 
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